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We study the localization and addressability of ultra cold atoms in a combined parabolic and 
periodic potential. Such a potential supports the existence of localized stationary states and we 
show that using a radio frequency field allows to selectively address the atoms in these states. This 
method is used to measure the energy and momentum distribution of the atoms in the localized 
states. We also discuss possible extensions of this scheme to address and manipulate particles in 
single lattice sites. 
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Periodic potentials have been used with great suc- 
cess in a series of experiments with ultra cold atoms 
01 For a trapped sample the periodic potential 
is usually accompanied with an additional parabolic con- 
finement 0? 3 • In most experiments involving trapped 
Bose-Einstein condensates only the properties of the 
ground state, which are not changed dramatically by the 
parabolic confinement, are important. However, for the 
excited states of the system the additional harmonic con- 
finement must be taken into account. This is especially 
true for atomic Fermi gases, where the Pauli principle en- 
forces a population of higher energetic states. The qual- 
itative different nature of the single particle states can 
be seen from recent experiments with ultra cold fermions 
which have found evidence for a localization of the parti- 
cles in such a potential In contrast to the Mott 
insulating phase |4j , this localization process is indepen- 
dent from the interaction and is a pure consequence of the 
potential shape. It is therefore important to understand 
its properties and physical consequences, nonetheless be- 
cause most of the theoretical work on quantum phase 
transitions involving fermions @, 0, OH UM are based on 
pure homogeneous systems. A combined periodic and 
parabolic potential is also interesting for possible applica- 
tions in quantum information and was recently proposed 
for the implementation of a qubit register for fermions 
[l^ . A possible addressability of individual atoms in sin- 
gle lattice sites is thereby an intriguing vision. 

In this work we study ultra cold atoms in a parabolic 
magnetic potential which is superimposed in its weak di- 
rection with a one dimensional optical lattice. The com- 
bined potential possesses two distinct classes of eigen- 
states, which — depending on their energy — either 
extend symmetrically around the trap center or are lo- 
calized on the sides of the potential. We use a radio 
frequency technique to address the particles in localized 
states and to measure the density of state along the lat- 
tice direction as well as the momentum distribution of 
the particles. Due to the localization of the particles, the 
radio frequency field addresses the particles in a denned 
spatial region and we discuss the possibility of extend- 



ing this scheme to manipulate individual atoms in single 
lattice sites. 

In the experiment we prepare an ultra cold cloud of 
87 Rb in the combined potential by forced evaporative 
cooling. The magnetic trapping potential has an axial 
and radial oscillation frequency of co a = 2ir x 16 Hz and 
uj v = 2tt x 197 Hz and the optical lattice (A = 830 nm) 
can be adjusted between < s < 10, where s mea- 
sures the potential height in units of the recoil energy 
E r = h 2 /2m\ 2 . The atoms are prepared in the spin- 
polarized | F = 2, mp = 2) state and the temperature 
is between 500 and 600 nK. Here we are only interested 
in the single particle behavior; in this frame bosons and 
fermions show the same features, as long as interparticle 
interactions can be neglected. 

To understand the properties of the combined potential 
we first solve the ID Schrodinger equation in the direction 
of the lattice, 
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where E n is the energy of the nth eigenstate. This Hamil- 
tonian has also recently been studied in tight binding 
approximation |l3|, UM- In Fig. la we show a density 
plot of the first 1000 eigenfunctions of QJ. Each line in 
Fig.^i corresponds to a density plot of the wave function. 
For low energies we find delocalized states that spread 
symmetrically around the potential minimum. Above a 
threshold energy, the wave functions of the eigenstates 
become localized on both sides of the potential. If we 
look at higher energies, a second group of eigenstates ap- 
pears, centered again around the trap minimum. It is 
straightforward to identify this shell-like structure with 
the well known band picture for a pure periodic poten- 
tial. This becomes particularly clear if one looks at the 
accessible energy values for a given position which cor- 
respond to the calculated bandwidth ^bw and band gap 
^gap of a pure sinusoidal potential, only shifted in energy 
by the local value of the parabolic potential. As a direct 
consequence, the system does not have an absolute but 
a spatially varying energy gap. The temperature of the 
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FIG. 1: (a) Spectrum of the Hamiltonian: representation of 
the ID spectrum of the single particle Schrodinger equation 
for a combined parabolic and periodic potential with s=3. 
Each line represents one eigenstate of the system, which is 
plotted as density profile in grayscale. The vertical position 
of the profile corresponds to the energy of the eigenstate. (b) 
Cloud of atoms without rf field (upper image) and with rf 
field (lower image) after 1.5 ms time of flight. 



cloud is chosen that the bandwidth of the first band is 
much smaller than the average energy of the atoms, thus 
providing a high population of the localized states. In the 
two radial directions, the particles occupy the harmonic 
oscillator states. 

We now introduce our experimental technique to pre- 
pare and address the particles in the localized states. Af- 
ter the end of the evaporation a radio frequency (rf) field 
is applied in order to induce spin flip transitions. Be- 
cause of the magnetic radial confinement this removes all 
atoms whose wavefunction has a spatial overlap with the 
magnetic field shell where the resonance condition hv = 
H&B(y)/2 is fulfilled. Periodically sweeping (1kHz rate) 
the radio frequency within an interval Av = v up — v\ ow we 
define a spatial region in which the atoms are removed 
from the potential. After 100 ms of rf field we image the 
atoms which are left in the potential. In Fig.^D we show 
an absorption image of the atomic cloud after 1.5 ms time 
of flight without and with rf field. In the latter case we re- 
main with two clouds, located at the edges of the original 
cloud. Atoms in these two clouds are trapped in localized 
states and the result directly shows, that the rf field is 
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FIG. 2: Energy distribution of localized states in a lattice 
with s=9. (a) Scan of the rf field through the cloud. The 
indicated frequencies are the upper frequency v up of the rf 
field, the sweep amplitude of the field is 3 kHz. (b) Atom 
number in the left cloud in dependence on is up . The straight 
line is a fit with an axial density of states that is proportional 
to E~ 1/2 (see text). 



capable to address the atoms in a defined spatial region. 
We have checked that (i) without the optical lattice all 
atoms are removed from the trapping potential and (ii) 
after switching off the rf field the atoms remain on their 
position in the trap. Only for small lattice heights we 
can observe a slow motion of the two clouds towards the 
trap center. This is due to presence of collisions which 
allow the bosons to hop between different localized states 
|6|. For a spin-polarized Fermi gas this effect is absent. 

In our recent work 0, Q we were able to detect the 
localization of the atoms by looking at the center of mass 
position of the whole cloud. Here, we can employ this 
kind of rf spectroscopy to look at the energy distribution 
of the atoms in the localized states. In Fig.J^t) we show 
a series of absorption images where we have scanned the 
rf field fixing the frequency interval Av = 3 kHz. In- 
creasing the rf frequency we start to remove atoms from 
the center of the trap. The hole in the spatial distribu- 
tion deepens, until the lower frequency bound is higher 
than the resonance frequency at the trap bottom: atoms 
in the center are no longer removed from the potential 
and we observe three clouds until for even higher fre- 
quencies the displaced peaks disappear and the cloud 
is again unaffected by the rf field. In Fig.[2h we show 
the number of atoms in the left cloud in dependence 
on upper frequency of the rf field. Due to the local- 
ization these atoms have an axial energy which is higher 
than E up = 2/i(^ up _ i , ), where v$ is the resonance fre- 
quency in the trap center. Thus, the number of atoms 
in the cloud is determined by the density of states in 
the axial direction: N ex n(I? a )p a (i? a )<iE a , where 
p a (i£ a ) is the density of states in the axial direction and 
n(E a ) = e - E */ k & T is the axial energy distribution. As 
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FIG. 3: Momentum distribution of the localized states. 

(a) Fourier transform of selected eigenstates of the spectrum 
shown in Fig.^i. The first three states are delocalized and 
lie within the first band, the last state is the lowest state of 
the second band. The states in between are localized states. 

(b) Absorption image of a cloud of atoms such as shown in 
Fig. 12b after 10 ms time of flight for s — 3 and s — 9. 



shown in Ref. 14] the density of states for energies larger 
than the bandwidth of the first band is predicted to be 
proportional to E~ x l 2 . Using this value we fit our data 
leaving the temperature and a constant of proportional- 
ity as free parameters 15]. The result which is shown in 
Fig.03 is consistent with our data. 

We now turn to the momentum distribution of the lo- 
calized states. For a potential depth of s = 3 we have 
calculated the Fourier transform of the wave functions of 
the eigenstates. Fig.|3Ji shows the momentum distribu- 
tion for selected eigenstates within the first and second 
band. The lowest eigenstate shows the well known peak 
distribution at multiples of twice the Bragg momentum. 
With increasing energy these peaks broaden and develop 
a substructure. The states of our interest are the local- 
ized states, whose distribution spreads over the first Bril- 
louin zone (±Pb) and is similar for all states, regardless 
the energy of the state. In the radial direction, the mo- 
mentum distribution is determined by the temperature 
of the cloud. Consequently, localized clouds as shown in 
Fig.0 are expected to exhibit an anisotropic expansion. 
In Fig.^b we show an absorption image of a localized 
cloud after 10 ms time of flight for s = 3 and 8 = 9. 
The measured aspect ratio of 2.5 for s = 3 is highly 
anisotropic and directly proves the nonclassical momen- 
tum distribution of the localized states. For s = 9 the 
cloud expands much faster in the direction along the lat- 
tice (horizontal direction). Indeed we calculate a 2 times 
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FIG. 4: (a) Tunnelling between the bands: density distribu- 
tion of the 594th eigenstate for a potential with s = 0.3. (b) 
Extension of the localized states in the first band in depen- 
dence on the lattice height for three different energies (/cb x 100 
nK, k B x 200 nK, and k B x 300 nK,). 



larger momentum distribution for s = 9 with respect to 
s = 3 which leads to a nearly isotropic expansion. 

We finish this work with discussing the localization 
process and its possible applications. For small lattice 
heights, the localization of the states is prevented by in- 
ter band transitions and as an example we show in Fig.^i 
an excited state for s = 0.3. The wave function exhibits 
substantial contributions from both bands and a particle 
in this state is no longer localized. On the other side the 
extension of the localized state shrinks with increasing 
lattice height (Fig.^b) and the smallest possible exten- 
sion is given by the ground state in each lattice site. For 
our parameters we find that for s = 30 the eigenstates 
are mainly located within a single lattice site. This re- 
sult is of particular interest because it shows that for 
similar experimental conditions like in the Mott insula- 
tor experiment Q a localization of the particles within 
one lattice site is possible without a repulsive interaction. 
Indeed, Viverit et al. jl^ have shown that if an atomic 
Fermi gas is loaded in a combined parabolic and periodic 
potential even an occupancy with exactly one atom per 
lattice site can be achieved. Another intriguing conse- 
quence of the localization is the addressability of single 
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lattice sites: the potential gradient discriminates the res- 
onance condition for an atomic transition in each lattice 
site if the transition depends on the external potential. 
In our setup the magnetic potential leads to a spatially 
varying Zeeman splitting within the \F — 2)-manyfold 
and thus a very weak radio frequency should allow - in 
principle - for the manipulation of the atoms within one 
lattice site. To get a reasonable discrimination and a 
sufficiently high Rabi frequency, the resonance condition 
between adjacent lattice sites should be shifted by some- 
thing about 10 kHz which would require a gradient of 
300 G/cm Thereby a linear potential is more fa- 

vorable than a parabolic one where the frequency shift 
is changing along the lattice. For well defined experi- 
mental conditions it would be also desirable to have an 
optical confinement in the radial direction because oth- 
erwise particles with lower axial but higher radial energy 
might also be resonant with the radio frequency. 

In conclusion we have proved that particles in a com- 
bined periodic and parabolic potential are trapped in lo- 
calized states. We have used a radio frequency field to 
induce spatially resolved transitions in order to remove 
the particles from the potential. Doing so it was possible 
to measure the axial density of states and the momentum 
distribution of the localized states. The experiment has 
implications in various directions. First it shows that 
an inhomogeneous periodic potential exhibit a qualita- 
tively different phenomenology compared to a homoge- 
neous system. Second the experiment directly evidences 
a new localization mechanism which is independent from 
the interaction. Especially for atomic Fermi gases such 
a potential can be interesting to create single lattice oc- 
cupancy [l2| . Third we show that in principle the parti- 
cles can be addressed spatially resolved. For future ex- 
periments it should be possible to extend the scheme to 
the manipulation of particles in single lattice sites which 
would constitute a major progress in "quantum engineer- 
ing" with ultra cold atoms. 
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